Does parton saturation at high density explain hadron muhiplicities at RHIC ? 
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We discuss the recent claim that hadron multiplicities measured at RHIC energies are directly 
described in terms of gluon degrees of freedom fixed from the initial conditions of central heavy ion 
collisions. The argument is based on the parton saturation scenario expected to be valid at high 
parton densities and on the assumption of conserved gluon number. Alternatively we conjecture 
that "bottom- up" equilibration before hadronization modifies this picture, due to nonconservation 
of the number of gluons. 



At RHIC and LHC energies in the central rapidity region of heavy ion collisions a high density of energy is deposited 
|. mainly in the form of gluons. In recent papers by Kharzeev et al. (hereafter referred to as KLN) it is assumed 

that the initial gluon density determines the number of produced (charged) particles, i.e. that there is a direct 
correspondence between the number of partons in the initial state and the number of particles in the final state. 
The argument is formulated in the framework of the saturation scenario [^|-^, which determines the initial gluon 
distribution inside the colliding nuclei, and where the characteristic momentum scale is the hard scale Qs » ^qcd- 
Based on this scenario the initial gluon multiplicity immediately after the high energy nuclear collisions can be 
calculated in the McLerran-Venugopalan model At tq ~ l/Q^ this initial gluon density is [Tl| ], 

C<\ ■ rihardiTo) = C — . 1 

^ , An-'Ncas [QsTo] 

' In this relation is used to obtain the hadron multiplicities in A-A collisions, which are then compared with 
, the corresponding RHIC data [^-16|. 



In this note we propose a possibly more realistic description, which departs from taking the initial condition as the 
I only one determining ingredient for the final hadronic state. This description is based on the "bottom-up" scenario 
[p7|-p^, which leads to thermalization of the gluons produced after the collision. In contrast to KLN, this scenario 
O I stresses the importance of branching processes of gluons which may allow for a short enough equilibration time, at 
^ least at high energies. As a consequence the number of gluons is increasing between initial and equilibration times. 
In this context let us note that already in po| a phenomenological factor Kinei > 1 is introduced accounting for gluon 
number changing processes which may occur at late times beyond when the classical approach is applicable. 

In the following we investigate in the "bottom-up" picture the resulting hadron multiplicities and compare with 
^ , the predictions of the KLN model. We discuss in some detail the conceptual differences between these two interesting 
• possibilities. 

The arguments of KLN go through a number of steps, which we first critically review: 

i) The number of gluons in the initial state, and at the time when gluons transform into hadrons, is assumed to be 
equal, which can be true when only 2 <-> 2 processes are taken into account. Only in this circumstance the measured 
hadron multiplicities are reflecting directly the initial conditions. 

ii) In (^ the parameter c is a constant linking the number of gluons in the nucleus wave function to the number of 
gluons which are freed during the collision. It is expected to be of 0(1) [Q. An (approximate) analytical calculation 
gives c = 2 ln2 ~ 1.39 |2l|. However, extracting c from the numerical simulation in |^,|oj leads to the estimate [] 
c ~ 0.5 for Q, = 1 GeV. 

iii) The hard saturation scale in the case of one nucleus has been determined to be , 

^ ^ ^\ "^(Qs) xG{x,Ql) Pnucleon{s) , (2) 



^This estimate is based on the relation c = 4:%'^ fN /[{N^ — 1) ^nQ^/AQ^j^], where /jv ~ 0.3 jiij 
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where sis the impact parameter, xG{x, Q"^) is the gluon structure function in the nucleon, and Pnucieon is the transverse 
density of nucleons in the nucleus. Notice that in the above equation an additional numerical multiplicative factor 
may come as a consequence of the inherent uncertainty in the precise determination of the saturation momentum 
outside of the McLerran-Venugopalan model We shall call this factor K and introduce it a little later on. KLN 
generalize (||) to the case of two colliding nuclei. Let us first write |^ 



where ppartis, b) is the density of participating nucleons in the transverse plane as a function of the impact parameter 
b of the A-A collision, and of s the transverse coordinate of the produced gluon. The factor i is required by the 
proper definition of Qs ; relative to one nucleus . Integrating the density with respect to s leads to 

d^S Ppart{s,b) = Npart[b) , (4) 

with Npart{b) the number of participants in the A-A collision at fixed b. 

iv) Focusing on the distribution of freed gluons dN/drj for 77 = at given b, we can write 

dN.-.. Nl^l f 1 

-(^)=c-rW / d's—Ql{s,b) 



d-q Att'^Nc J a. 



=,2^ Npart{b) 



cxG(x,Qf)i:i^. (5) 



The corresponding charged hadron multiplicity for the most central collisions is then obtained as 



Npart dri 3 



In 



^QCD 



(6) 



using the gluon structure function (||) given below. As reference energy we take ^/s = 130 GeV. 

One should remark that in (||) and shows up as an effective average over the variable s in (^. As a first 

approximation the following is used: 

Ql(&) = K a.m xG{x, Ql) . (7) 

v) Information on the gluon structure function is necessary. In the small x regime, the main feature is that it 
increases with at fixed x. Following KLN it is reasonable to take 



xG(a:,Q2) ^0.5 ln(-^), (8) 

^QCD 



with Aqcd = 200 MeV, such that xG{x, Q^) ~ 2 at ^ 2 GeV^ (at x = 0.02). 
vi) The strong coupling constant is 



"^(Q') ^ , . / . ^ (9) 
f3o ln(^ 



QCD 



with /3o = (11 — 2ny/3)/47r for = 3; we take Uf = ?>. 

vii) Using (|^) together with the above choices for as and xG{x,Q^), one finds for central Au — Au collisions at 
RHIC (at = 130 GeV), taking for the moment K = I: 

Q^^{b = 0) ~ 0.63 GeV^, using Ppart{b ^ 0) = 3.06 fm-^, 



^Notice in (^, consistently with (P), the factor i compared to eq. (14) in [Q]. This has the consequence of lower values of 
appearing in our discussion of the KLN approach. 
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as quoted in Table 2 by For larger values of 6, becomes even smaller, e.g. for 6 = 10 /m, Qs — 0-32 GeV"^ ! 
This casts some doubts on the applicability of this model, even for central collisions at RHIC energies. 

As already noticed, the missing factor ^ in KLN and a larger value used for a^, namely as{Q^) = 0.6 at Q'^ = 2 GeV^, 
allows them to get larger values of Q^, i.e. > 2 GeV^ for the central collisions. 

Taking in the following an optimistic point of view, we use with a multiplicative K factor of order 0(1), explicitely 
K ~ 1.6, such that for 6 = 0, Q^ib = 0) = 1 GeV^. Our use of if is to a large extent cosmetic. This factor only 
appears in calculating Ql , but does not change (||) or (|l8|) . Such a factor will affect the average transverse momentum 
per produced gluon but not the total number of produced gluons. The difference between Qs{b = 0) = 1 GeV and 
Qs{b = 0) = 0.8 GeV, corresponding to K = 1.6 and K = 1, respectively, has little effect on the equilibration 
temperature Teq and time Teq quoted later in this note. 

viii) Finally, there is an equality assumed between the number of partons in the final state and the number of 
observed hadrons ("parton-hadron duahty" p^). 

Let us now turn to the "bottom-up" scenario Jl7|-p^, which in contrast to the KLN ||l|-|| prescription, does not 
relate the multiplicities to the initial condition only, but also to the way gluons are thcrmalizcd. 

In the framework of perturbative QCD the time evolution of the gluonic system, when described by a non-linear 
Boltzmann equation based on 2 2 processes, a relatively long time (~ Qj^ exp(Q;s ^^^)) is required for the 
approach to kinetic equilibration pl],E4|. Fast and efficient thermalization occurs when inelastic processes, namely 



gluon splittings 2 3 are taken into account, and kinetic equilibration occurs much faster at times ~ as ^^^^Q^ ^ 
jiot . In this "bottom- up" scenario the time evolution of the system proceeds through several regimes (Fig. |^), 



with QsT ^^"^ -^"^ ^'^^^ 



The difference between the "bottom- up" and the KLN picture is schematically illustrated in this Fig. |l]. Under 
the KLN assumption the hard gluons (on the momentum scale Qs) are conserved in number, i.e. UhardT — const, 
and they finally hadronize, after passing through a hydrodynamical stage. The "bottom- up" scenario is characterized 
by the fact that hard gluons are degrading, soft ones are formed and start to dominate the system. As a result the 
interactions of gluons in this kinetic scenario modify strongly the initial gluon spectrum. The gluons are redistributed 
and thermalizing, such that a quark gluon plasma is formed. The number of gluons, together with the entropy [ p^ , 
is increasing with proper time r, such that the ratio is, 

R = [nsoMT){QsT)]\r^, I {nharArmsT)\U ~ "J^/S » 1 . (10) 

The following processes are expected to be present (at RHIC and higher energies in the central region of pseudorapidity 
7y < 1): at QsT > 1 saturated hard gluons — > elastic scatterings and branching/production of soft gluons — > at 

QsT\eq ~ Ois ^^^^ thcrmalizatiou of these soft gluons with temperature Teq ~ ai^^Qs — > hydrodynamic expansion 
— > hadronization (Fig. 

In more detail the parametrically estimated time scales are determined as follows [p7|-p^ . 

Early times 1 <C QsT ^ as • 

At the earliest time, t ~ Q^^, gluons, i.e. hard gluons, have typical large transverse momentum of order Qs 
and occupation number of order 1/as. Later on gluons with smaller momenta, but still larger than Aqcd will be 
produced (nevertheless, they are denoted as soft gluons). The density of hard gluons Uhard decreases with time due 
to the one-dimensional expansion. Gluons interact by elastic scatterings at small angle, with exchange momentum 
<C Qs- The typical occupation number is large until QsT ~ as ^^'^ , when it becomes of 0{1). 

This regime is the transition region from the non-linear classical gluon field to the one where the transport description 
by Boltzmann equations should become applicable. 

Times as"^/^ < QsT < a^^^^ 

1 /2 

Inelastic scatterings produce (soft) gluons Usoft with characteristic momentum estimated to be of order as Qs > 
Aqcd (Fig. P, namely via hard+ hard — > hard+ hard + soft. The number of these soft gluons becomes comparable 

to that of hard ones at QsT ~ as ^^'^ , namely risoft ~ rihard- 
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FIG. 1. Characteristic momentum scales for the "bottom-up" scenario. 



Times QsT 3> 

After QgT ~ as most gluons are soft, Ugoft 3> rihard', they achieve thermal equihbration amongst themselves. 
Although the whole system is still not in thermal equilibrium, the soft part is characterized by the temperature T, 
nsoft ^ ■ The few hard gluons collide with the soft ones of the thermal bath and constantly loose energy to the 
latter. A hard gluon emits one with a softer energy, which splits into gluons with comparable momenta. The products 
of this branching quickly cascade further, giving all their energy to the thermal bath. This increases the temperature 
in the thermal bath found to be. 



T 



CT alQlr , 



(11) 



i.e. it increases linearly with time, even when the system is expanding, due to the hard gluons which serve as an 
energy source, ct is a numerical constant to be discussed later. One can verify a posteriori that soft gluons are indeed 
equilibrated due to many interactions: the size of the expanding system of 0{t) is indeed much larger than the mean 
free path of the soft gluons, i.e. 



(12) 



where the cross section is estimated as ct ~ a^/T'^. 

In the following, in order to provide predictions for particle multiplicities in this scenario we go as much as possible 
beyond the parametric estimates, ending up with consistency arguments which specify the allowed range of parameters. 

First we start by considering the constant ct which may be written in terms of the parameter c as [ p7| . 



CT 



15 



0.16 c. 



The linear growth of T as shown in Fig. |l| terminates, when the hard gluons loose all of their energy, i.e. when 



eq 



-13/5. 



(13) 



(14) 



where the parameter Ceq is unknown for the moment, although in principle it could be calculated in the "bottom- up" 
framework. The temperature achieves a maximal value, i.e. the equilibration temperature, which is expressed by 
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Teg = 0.16 c ce, a^/'(Q^) Q 



(15) 



Subsequently the temperature decreases as r ^/"^ p5| , such that rrisoftiT) = const. 

We now derive the charged hadron muhiphcity from the number density of the equihbrated soft gluons, 



2 i\C(3)y3 



eq ' 



which is much larger than Uhard in (|l|) used by KLN. One finds for the ratio R defined in (JK 

For the charged hadron multiplicity the result for = 130 GeV and the most central collisions, is 

, 2 dA^c/i, Rc 



part 



drj 



In —77 
3 A2 



0.04 c 



,3 „4 



In 



^QCD 



7/5 



which replaces 



(16) 



(17) 



(18) 



In order to compare the predicted charged hadron multiplicity, (|^) and (18), with RHIC data we use as a reference 
the result by the PHOBOS Collaboration |]l^, namely 3.24±0.1(siai)±0.25(si/si), which is in good agreement within 
errors with the experimental measurements of the other collaborations at RHIC. This is best seen by comparing the 
values for dNch/drj at midrapidity at = 130 GeV: dNch/dri = 555±12{stat)±35{syst) 609 ±l{stat)± 37 (syst) 
0, 5i9±l{stat)±35{syst) 567 ± l{stat) ± 38{syst) @, respectively. 
As reference value we take 



2 dN,h^ 



part 



drj 



3.24. 



(19) 



First we note that in the KLN approach agreement with experimental data can only be achieved with the value 
c ~ 3, which is different from the current numerical estimate quoted before, c ~ 0.5. 

Turning to the "bottom- up" approach, the experimental value (|9|) meets the theoretical expectation (18) for 
i? c ~ 3, or equivalently 



2.0 



(20) 



This relation is to be confronted with the consistency requirement for the "bottom- up" scenario , that the ratio R 
( |l7| ) be larger than 2, implying 



ctq > 10, i.e. Ceq > 



(21) 



It is not too difficult to see that (EG, pl^) constrain the two parameters: 



c < 1.5 and Ceq > 1.5. 



(22) 



This in turn allows us to infer the actual properties of the medium, especially to answer tentatively the question of 
the formation of the equilibrated plasma. We do this via discussing the temperature T^q and the equilibration time 
Teq. In order that the quark gluon plasma is produced, T^q should be bigger than the phase transition temperature 
Tdeconf, which is of order Tdeconf — 173 ± 8 MeV and 154 ± 8 MeV for 2 and 3 flavour QCD, respectively i.e. 



^eq 2^ ^deconf ■ 



(23) 
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This constrains 



c Cen > 1.3, or c > 0.2 . 



(24) 



Finally, we may correlatively discuss Tgq . Under the condition ( p2| ) one finds (for central collisions) 

Teg > 2.6 fm^ 1/2 Rau, 



(25) 



which is much bigger than the current estimate of ~ 0.7 fm |^ Within the uncertainties inherent to these estimates 
the formation of the equilibrated plasma may indeed be realized. 

One may consider the energy dependence at RHIC energies of the charged multiplicities. Following |^ it is controlled 
by the energy dependence of the saturation scale, i.e. Q'^{s)/Q'^{so) = (s/sq)^^'^ with A = 0.25 |28|. As discussed we 
choose as the reference value Qs(so) = 1 GeV^ at = 130 GeV. For the most central collisions the expression for 
the energy dependence of the pseudorapidity density of charged particles at midrapidity reads in the "bottom- up" 
scenario, 



2 dNch 



) ~ 0.64 



"part drj 

Based on ( p6| ) the multiplicity at -^s — 200 GeV is obtained: 

, 2 dNch 



In 



^^QCD 



7/5 



part 



drj 



= 3. 



compared to the value given by the PHOBOS Collaboration [|r 

, 2 dNeh 



part 



drj 



= 3.78 ±0.25 (syst). 



(26) 



(27) 



(28) 



Finally we consider the centrality dependence, i.e. the dependence of dNch/di] as a function of Npart- In Fig. 
the comparison of the "bottom- up" expectation as derived from ( p6| ) with data from the PHOBOS Collaboration |T 
at ^/s = 130 GeV and at = 200 GeV at RHIC is shown. For this comparison Q^(6) is calculated from 



as a 



function of Npart, using Ppart{b) given in Table 2 of ||l|, together with the scaling relation for the energy dependence 
of Qs quoted above pM]. One has to note, that for Npart < 100 the values of Ql are becoming smaller than 0.6 GeV^. 



The structure of the shape of dNch/drj as a function of Npart seen in the data (Fig. 
(^, to details of xG{x, Ql) at small function of Q^. 



3) could be attributed, cf. with 



^This short equilibration time is required in order to describe the hadron spectra for different particles as measured at RHIC 
by hydrodynamic calculations, assuming the existence of the quark gluon plasma . 
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FIG. 2. The scaled pseudorapidity density for |r;j < 1 as a function of Npart at ^/s — 130 GeV and 200 GeV. Data are from 
the PHOBOS Collaboration, and curves from the "bottom-up" scenario, as described in the text. 



In summary, the description provided by the "bottom-up" scenario is in agreement with RHIC data, provided the 
parameters c and Cg^, which are not determined in the picture, lie in a given, limited range. In particular for c ~ 1 and 
Ceq — 2 the picture looks both reasonable and attractive. For these values the results (for the most central collisions) 
for R,Teq and Teq are 3, 230 MeV and 3.6 /m, respectively, for — 1 GeV; for Qs — 0.8 GeV, i.e. for K ~ 1, the 
values are i? = 3, T^q = 210 MeV and T^q = 3.2 fm. However, the picture probably does not make much sense for 
a small value of c ~ 1/2, which is currently favoured by numerical calculations of classical field equations, because 
the ratio R turns out to become i? ~ 6, which is too much inelasticity when Qs — I GeV. This discrepancy deserves 
further examination. 

On the other hand, in order to accomodate the present RHIC data, the KLN description requires c ~ 3. For the 
moment, due to a number of ambiguities, it is not yet clear whether this large value is compatible with the various 
constraints of the saturation picture. 

The "bottom-up" scenario is likely to provide a more convincing agreement with data. However, this analysis will 
have to be supplemented by further ingredients, e.g. the calculation of c^q, a more precise estimate of c, before one 
can finally claim a significant agreement with data. 
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